Unmasking the nuclear matter equation of state 
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Accurately calibrated (or "best fit") relativistic mean-field models are used to compute the dis- 
tribution of isoscalar monopole strength in ®'^Zr and ^°*Pb, and the isovector dipole strength in 
^''^Pb using a continuum random-phase-approximation approach. It is shown that the distribution 
of isoscalar monopole strength in '^"^Pb — but not in ^"Zr — is sensitive to the density dependence of 
the symmetry energy. This sensitivity hinders the extraction of the compression modulus of symmet- 
ric nuclear matter from the isoscalar giant monopole resonance (ISGMR) in 2°**Pb. Thus, one relies 
on ^°Zr, a nucleus with both a small neutron-proton asymmetry and a well developed ISGMR peak, 
to constrain the compression modulus of symmetric nuclear matter to the range K= (248 ±6) MeV. 
In turn, the sensitivity of the ISGMR in ^"^^Pb to the density dependence of the symmetry energy 
is used to constrain its neutron skin to the range 7?„ — 7?p<0.22 fm. The impact of this result on 
the enhanced cooling of neutron stars is briefly addressed. 

PACS numbers: 24.10.Jv, 21.10.Re, 21.60.Jz 



Constraining the equation of state (EoS) of neutron- 
rich matter remains a fundamental problem in nuclear 
physics and astrophysics. The stability of neutron-rich 
nuclei [ij , the dynamics of heavy- ion collisions 0, 13 , the 
structure of neutron stars Q , and the simulation of core- 
collapse supernova H, 0|i ^ depend sensitively on the 
EoS. Unfortunately, our window to the EoS is limited 
by terrestrial experiments that have, until now, probed 
stable nucleonic matter at (or close to) nuclear-matter 
saturation density. Fortunately, dramatic improvements 
are unfolding on several fronts. First, the commissioning 
of new radioactive-beam facilities all over the world will 
probe the EoS at large neutron-proton asymmetries. By 
defining the limits of nuclear existence, these exotic nuclei 
will constrain the EoS of neutron-rich matter at (and be- 
low) normal nuclear densities. Second, space-based tele- 
scopes have started to place important constrains on the 
high-density component of the EoS New telescopes 

operating at a variety of wavelengths are turning neutron 
stars from theoretical curiosities into powerful diagnostic 
tools. 

The nuclear matter equation of state is conveniently 
parametrized in terms of the energy of symmetric nuclear 
matter {B/A) and the symmetry energy {S/A) in the 
following form: 
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at a density of po — 0.15 fm~^ (fc^ ~ 1.3 fm~^) and a 
binding-energy per nucleon of eq — — 16 MeV. It should be 
noted that one of the main virtues of the above Taylor- 
series expansion around saturation density [Eq. (^] is 
that the linear term in ^ for symmetric nuclear matter 
{i.e., the pressure) automatically vanishes. Yet no such 
special saturation point exists in the case of the symmetry 
energy. Indeed, the symmetry energy at saturation den- 
sity is not well known. Rather, it is the symmetry energy 
at the lower density of po - 0.10 fm"^ (fc^ ~ 1.15 fm"^) 
that seems to be accurately constrained (to within 1 
MeV) by available ground-state observables 0, -"-t 
should be emphasized that present-day experiments can 
fix only one isovector quantity. If one insists — and one 
should not — on constraining the parameters of the sym- 
metry energy at saturation density, then one would find 
a str ong correlation among its parameters (J, L, ifgymj 
. . .) |l3- For example, relativistic models consistently 
predict larger values for both the symmetry-energy coef- 
ficient J and the slope L at saturation density relative 
to nonrclativistic Skyrme models. This must be so if all 
models are to reproduce the value of the symmetry en- 
ergy at the lower Fermi momentum of kPp. Thus, in the 
present contribution we adopt the following convention: 
the symmetry energy is expanded around fc^ ~ 1.15 fm^^ 
and the value of the symmetry energy at fc^ is fixed at 
J = 25.67 MeV. That is, 



SIA{kY) = J - 
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Here the deviation from the equilibrium Fermi momen- 
tum is denoted by ^= (fcp - fcF)/fcF' tl^e neutron-proton 1 e M 7o^/7o M 4- 4-1, + I, f 4.1, ia-i^ " 
asymmetry by 6 ^ {N-Z)/A, and the various coefficients ^ = (^/ " ^F)/fcF- Note that henceforth tilde 



{K , J, L, Ksyra) parametrize the density dependence of 
the EoS around saturation density. 

Seven decades of nuclear physics have placed impor- 
tant constraints on the nuclear matter equation of state. 
Indeed, the energy systematics of medium to heavy nu- 
clei, when combined with accurately calibrated models, 
place the saturation point of symmetric nuclear matter 



quantities refer to parameters of the symmetry energy at 
fc^~1.15 fm"^ 

Having established that existing ground-state observ- 
ables accurately determine the binding energy per nu- 
cleon eg at k^p and the symmetry-energy coefficient J at 
k'p , how can one constrain any further the density depen- 
dence of the equation of state? In the case of symmetric 



2 



nuclear matter, the dynamics of small density fluctua- 
tions around the saturation point is controlled by the 
compression modulus K . The isoscalar giant monopole 
resonance (ISGMR) in heavy nuclei has long been re- 
garded as the optimal observable from which to deter- 
mine the compression modulus llj. This is especially 
true now that the breathing mode has been measured on 
a variety of nuclei with unprecedented accuracy jl2j . In 
contrast, the density dependence of the symmetry en- 
ergy is poorly constrained. Indeed, one may fit a va- 
riety of ground-state observables (such as charge densi- 
ties, binding energies, and single-particle spectra) using 
accurately-calibrated models that, nevertheless, predict a 
wide range of values for the neutron skin of ^°^Pb • As 
the neutron skin of a heavy nucleus is strongly correlated 
to the slope of the symmetry energy ^9., lOj. measuring the 
skin thickness of a single heavy nucleus will constrain the 
density dependence of the symmetry energy. The Parity 
Radius Experiment (PREX) at the Jefferson Laboratory 
aims to measure the neutron radius of ^"^Pb accurately 
(to within 0.05 fm) and model independently via parity- 
violating electron scattering |Q, ^|. This experiment 
should provide a unique observational constraint on the 
density dependence of the symmetry energy. 

While the above arguments suggest a clear path to- 
ward constraining the density dependence of the EoS, 
theoretical uncertainties have clouded these issues. First 
and foremost is the apparent discrepancy between non- 
relativistic and relativistic predictions for the value of 
the compressional modulus of symmetric nuclear matter 
required to reproduce the ISGMR in ^°*Pb. While non- 
relativistic models predict 220-235 MeV [Tlll7lll8l|. 
relativistic models arg ue for a significantly larger value 
if ~ 250-270 MeV ElSim. Further, relativistic mod- 
els systematically predict larger values for the neutron 
skin of 208p]-, i-giative to nonrelativistic Skyrme models. 
One of the goals of this contribution is to show that these 
two points are related. Indeed, the aim of this contribu- 
tion is twofold. First, to vindicate — through the exclusive 
use of accurately-calibrated models — our previous asser- 
tion that the distribution of ISGMR in heavy nuclei, and 
therefore the inferred value oi K , is sensitive to the den- 
sity dependence of the symmetry energy [23 |. Second, 
to rely on existing data on the isoscalar giant-monopole 
resonance in ^°Zr and "^^^Vh ^"^^ the isovec- 

tor giant-dipole resonance in 208p|-| j2^^ limits — 

simultaneously — on the compression modulus of symmet- 
ric nuclear matter and on the neutron skin of ^°*Pb. 
Note that since first proposed , other groups have ad- 
dressed the possible impact of the density dependence of 
the symmetry energy on the ISGMR in ^ospb [Uiim. 

The starting point for the calculations is an interacting 
Lagrangian density of the following form ^13, 26]: 



where $ = gs<^, = g^V^j., and = gph^. The La- 
grangian density includes Yukawa couplings of the nu- 
cleon field to a scalar (</>) and to three vector (V^^ , b'^, 
and A^) fields. In addition to the Yukawa couplings, 
the Lagrangian is supplemented by three nonlinear me- 
son interactions. The inclusion of scalar- meson interac- 
tions (via K and A) is used to soften the equation of state 
of symmetric nuclear matter, while the mixed isoscalar- 
isovector coupling (Av) modifies the density dependence 
of the symmetry energy — without affecting well known 
ground-state properties. Note that this last term was 
absent from Ref. so the softening of the symmetry 
energy had to be done artificially. This drawback has 
now been corrected. 

The relativistic mean-field models employed here 
are motivated by the enormous success of the NL3 
parametrization [13. Indeed, the NL3000 set used here 
(having Av = 0) is practically identical to the original 
NL3 model. The other sets are obtained by adding an 
isoscalar-isovector coupling Av 7^ 0, while at the same 
time re-adjusting the strength of the NNp coupling con- 
stant (gp) to maintain the symmetry-energy coefficient 

fixed at J = 25.67 MeV (see discussion above). The 
aim of this added coupling is to change the neutron den- 
sity of heavy nuclei, while leaving intact ground-state ob- 
servables that are well constrained experimentally. One 
should stress that the addition of Av has no impact on 
the properties of symmetric nuclear matter, so the sat- 
uration properties remain unchanged. In summary, all 
the models used in this contribution share the following 
properties with the original NL3 model of Ref. : for 
symmetric nuclear matter, a Fermi momentum at satu- 
ration of kp = 1.30 fm~^ with a binding-energy per nu- 
cleon of £0 = —16.24 MeV, and a compression modulus 
of K — 271 MeV. For the symmetry energy, a symmetry- 
energy coefficient of J = 25.67 MeV at a Fermi momen- 
tum of A;^ = 1.15 fm^^. 




.930- (5v^^^ + f r • b^ + |(l + r3)A^ 



- + Av {W^Wn {Bp ■ B^) , (3) 



FIG. 1: Distribution of isoscalar monopole strength in ^"Zr for 
the small momentum transfer of g = 45.5 MeV. The response 
includes a small artificial width of 0.5 MeV. 
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While the success of the NL3 interaction in reproduc- 
ing ground-state properties (such as binding energies, 
charge radii, energy separations, etc.) for a variety of nu- 
clei all throughout the periodic table is well documented, 
we display in Table ^ ground-state properties for only 
the two nuclei of relevance to this contribution, namely, 
^°Zr and ^°*Pb. However, even these accurately cali- 
brated models predict a wide range of values for the neu- 
tron skin of ^°^Pb, confirming that the neutron-skin of 
a heavy nucleus is not tightly constrained by known nu- 
clear observables. The fifth column in the table displays 
the compression modulus of asymmetric nuclear matter 
with a neutron-proton asymmetry corresponding to ^^Zr 
(6 = 0.111) and ^°^Ph {b = 0.212). It is this quantity— not 
the compression modulus of symmetric nuclear matter — 
that is constrained by the breathing mode of nuclei. This 
simple fact makes the connection between the measured 
ISGMR and the compression modulus of symmetric nu- 
clear matter sensitive to the density dependence of the 
symmetry energy. Recall that the compression modulus 
of symmetric nuclear matter was fixed in all models at 
K = 271 MeV, yet for {b = 0.212) asymmetric nuclear 
matter the compression modulus ranges from 243 MeV 
(for the stiffest symmetry energy) all the way up to 
260 MeV (for the softest symmetry energy). Finally, the 
last column in the table shows peak and centroid ener- 
gies for the ISGMR in ^'^Zr and ^"^Pb computed in a 
relativistic random-phase-approximation (RPA); Figs.Q] 
and|21display the corresponding distribution of strength. 
As expected, there is a strong correlation between the 
centroid energies and the compression modulus of asym- 
metric nuclear matter. One should note in passing that 
the continuum RPA formalism employed here, but re- 
ported elsewhere [23, , respects important symmetries 
of nature, such as translational invariance (in the form of 
Thouless' theorem |2^|33) and the conservation of the 
vector current. 




CO (MeV) 

FIG. 2: Distribution of isoscalar monopole strength in ^°*Pb 
for the small momentum transfer of g = 45.5 MeV. The re- 
sponse includes a small artificial width of 0.5 MeV. 



The great advantage of a nucleus such as Zr is that it 
has both a well developed isoscalar-monopole peak and 
a small neutron-proton asymmetry (6 = 0.111). The lat- 
ter manifests itself into the near collapse of all curves in 
Fig. ^ into a single one, so that the former may directly 
constrain the compression modulus of symmetric nuclear 
matter. In contrast to ^°Zr, the distribution of ISGMR 
strength in ^^^Pb is sensitive to the density dependence 
of the symmetry energy. While this sensitivity should be 
sufficient to constrain the density dependence of the sym- 
metry energy, one could do even better. Indeed, one may 
constrain the density dependence of the symmetry energy 
by demanding that both the ISGMR and the isovector 
giant-dipole resonance (IVGDR) in 208p|-, simultane- 
ously reproduced. The distribution of isovector dipole 
strength in ■^°*Pb is displayed in Fig.|31 We note that the 
isovector-dipole response gets hardened as the symmetry 
energy is softened. As all models share the same value of 
the symmetry-energy coefficient at A;^ = 1.15 fm^^, the 
hardening of the response follows as a result of the sym- 
metry energy being higher at the (low) densities relevant 
to the isovector-dipole mode i21il . 
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FIG. 3: Distribution of isovector dipole strength in Pb for 
the small momentum transfer of g = 45.5 MeV. The response 
includes a small artificial width of 0.5 MeV. 

To constrain simultaneously the compression modu- 
lus of symmetric nuclear matter and the neutron ra- 
dius of ^'^^Pb, one starts by noticing that the theoreti- 
cal centroid energy of the ISGMR in ^°Zr overestimates 
the experimental value by about 1 MeV. Although a 
proper adjustment of K should be done through a re- 
calibration of parameters, a simple, yet accurate esti- 
mate may be obtained via the following scaling relation: 
^'ISGMR oc ^/K Using this relation and accounting 
for experimental uncertainties, an adjustment of about 
20 MeV in K is required to reproduce the ISGMR in 
9°Zr. That is, if = 271 MeV ^ = (248 ± 6) MeV. This 
adjustment in K induces a corresponding correction in 
the calculated values of the ISGMR in ^ospb. The cen- 
troid energies after correction, together with the peak 
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TABLE I: Binding energy per nucleon, root-mean-square charge radius, neutron-minus-proton root-mean-square radius, com- 
pression modulus for asymmetric (6 = 0.111 and) nuclear matter, and peak and centroid GMR energies for ^"Zr in the various 
models discussed in the text. The binding energy includes a center-of-mass correction of —0.08 MeV/nucleon, while the centroid 
energy (mi /mo) was computed by generating the distribution of strength in the range 10 < lj < 26 MeV. The second set of 
numbers in the table are for ^°®Pb (6 = 0.212) with a center-of- mass correction of —0.02 MeV/nucleon and a centroid energy 
extracted from a distribution of strength generated in the range 8<uj<23 MeV. 
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energies of the IVGDR in ^°^Pb, are displayed in Fig. El 
alongside the experimental values 0, [23 . The numbers 
in parenthesis indicate the predicted values for the neu- 
tron skin in ^°*Pb. The figure suggests that models with 
neutron skins in ^°^Pb larger than Rn — Rp — 0.22 fm 
may have an unrealistically stifT symmetry energy. 
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FIG. 4: Comparison between theoretical and experimental 
ISGMR centroid and IVGDR peak energies for ^°®Pb. Quan- 
tities in parenthesis represent the predictions for the neutron- 
skin of 2°**Pb in the various models discussed in the text. 



In summary, relativistic mean-field models have 
been used to compute the distribution of isoscalar- 



monopole strength in ^"Zr and ^°*Pb, and of isovector- 
dipole strength in ^"^Pb using a continuum RPA ap- 
proach. It was demonstrated — using exclusively ac- 
curately calibrated models — that the distribution of 
isoscalar monopole strength in ^°^Pb is sensitive to the 
density dependence of the symmetry energy. Further, 
existing experimental data were used to set limits on 
both the compression modulus of symmetric nuclear mat- 
ter and on the neutron skin of ^°^Pb. It appears that 
medium-mass nuclei, having a well-developed ISGMR 
peak and a small neutron-proton asymmetry (such as 
^"Zr but not ^°*Pb) allow for the best determination of 
the compression modulus of symmetric nuclear matter. 
In turn, the sensitivity of the ISGMR and the IVGDR in 
^•^^Pb to the density dependence of the symmetry energy 
may be used to impose constraints on the neutron skin 
of ^°^Pb. From the present analysis, a compression mod- 
ulus of K= (248 ± 6) MeV and a neutron skin in ^ospb 
of Rn—Rp ^ 0.22 fm were obtained. These values appear 
closer to those predicted in nonrelativistic studies. 



We conclude with a comment on the impact of these 
results on the cooling of neutron stars. In earlier publi- 
cations we have demonstrated how improved values for 
neutron radii could have a widespread imp act on the 
structure and dynamics of neutron stars [la. l2a . l3lL l33 | . 
In particular, we suggested that the enhanced cooling of 
the neutron star in 3C58 33] may be due to the conven- 
tional URCA process — provided the symmetry energy is 
stiff enough to generate a neutron skin in ^°*Pb larger 
than Rn — Rp > 0.24 fm. In view of our present findings, 
this now seems unlikely. Thus, the possibility that 3C58 
harbors an exotic star, such as a quark star, looms large. 
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